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Abstract 

We show the necessity of two-nucleon axial currents and associated pion 
emission/absorption operators for the partial conservation of the axial cur- 
rent (PCAC) nuclear matrix elements with arbitrary nuclear dynamics de- 
scribed by a nonrelativistic Schrodinger equation. As examples we construct 
such nonrelativistic axial two-body currents in the linear- and the heterotic 
{gA = 1-26) sigma models, with an optional isoscalar vector (u) meson ex- 
change. The nuclear axial current matrix elements obey PCAC only if the 
nuclear wave functions used in the calculation are solutions to the Schrodinger 
equation with the static one-meson-exchange potential constructed in the re- 
spective (sigma) model. The same holds true for the nuclear pion production 
amplitude, since it is proportional to the divergence of the axial current ma- 
trix element, by virtue of PCAC. Thus we found a new consistency condition 
between the pion creation/absorption operator and the nuclear Hamiltonian. 
We present examples drawn from our models and discuss the implication of 
our results for one-pion- two-nucleon processes. 

PACS numbers: 25.30.-c, 25.80.Hp, 25.10,+s 
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I. INTRODUCTION 



In an earlier publication pj a systematic study of axial current (partial) conservation 
was begun in the Bethe-Salpeter (BS) approach to nuclear bound states. There it was 
found that partial conservation of axial current (PCAC) puts constraints not only on the 
form of the axial current operator, but also on the nuclear wave function, by way of fixing 
the "potential" entering the nuclear BS equation. Since the pion production/absorption 
amplitude is an integral part of the conserved axial current nuclear matrix element, the 
same constraints are imposed on it as well. It stands to reason that the same kind of 
constraint will carry over into the nonrelativistic (NR) formalism. [] In this paper we extend 
the work in Ref. |l[ to the nonrelativistic description of the nucleus. We do not make a 
direct nonrelativistic reduction of the BS equations and amplitudes from Ref. |IJ, for they 
are well known not to have a good NR limit, but rather use this reference as a guide to 
developing the corresponding result within the Schrodinger equation approach. We shall 
show that even at this nonrelativistic level there are differences between the various versions 
of the a model used in Ref. 0]. 

The aforementioned PCAC consistency condition between the (elementary) pion pro- 
duction mechanism and the NN potential determining the nuclear wave equation is a novel 
feature promising to introduce a higher level of logical coherence into present day calculations 
of pion production on nuclear targets. As an example one may consider the pp — > 7c°pp and 
pp — > ir + d reactions. There, in particular, the 7r-production mechanism (operator) has often 
been considered without any reference to the nuclear wave functions 0|§. Several models 
involving scalar a-meson-exchange pion production "current", besides the 7r-exchange one, 
have been proposed [||,[5] . We shall show here exactly which ex- meson-exchange 7r-production 
operator must be included and how when the two-nucleon potential contains a one-sigma- 
exchange term, and vice versa. It turns out that this 7r-production operator is rather sensitive 
to the details of the sigma model used. We also analyze the tu-meson-exchange 7r-production 
operator required by the presence of the cj-meson-exchange potential. 

This paper falls into five sections. After the Introduction, in Section II, we define the 
context of our analysis and prove our most general result: the necessity of two-body axial 
currents for interacting two- (or more) particle systems obeying nonrelativistic quantum 
mechanics. In section III we construct two-nucleon axial currents that respect PCAC at the 
level of nuclear matrix elements, starting from three underlying relativistic chirally symmet- 
ric meson-nucleon model Lagrangians. In Section IV we present the three corresponding 
sets of two-nucleon 7r-production operators as an example of the main proposition and we 
discuss the results. In Section V we summarize and draw the conclusions. 



The idea that PCAC determines the pion-nuclear production operator has been around at least 
since the work of Blin-Stoyle and Tint . But the notion that the nuclear wave functions entering 
the same pion production amplitude are also constrained by the PCAC appears to be new. 
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II. PARTIAL CONSERVATION OF NUCLEAR AXIAL CURRENT 



The notion of PCAC is both historically and conceptually the foundation of chiral sym- 
metry in hadronic interactions. Modern implementations of this symmetry, such as the 
nuclear chiral perturbation theory (NxPT), however, do not emphasize that point. In the 
following we shall try and present another viewpoint of PCAC that will bear significance to 
nuclear physics in general and pion-nuclear processes in particular. In short, PCAC states 
that the hadronic axial current J^ 5 must satisfy the following continuity equation 

d'JU = -Kn a , (i) 

or equivalently 

V.J^(R) + ^^ = -/.m^(R), (2) 

where IP is the (canonical) pion field operator. In the quantum mechanical framework 
this can be written as an equation relating the divergence of the three-current and the 
commutator of the Hamiltonian and the axial charge density: 

V • Jg(R) +i[H, p a 5 (R)\ = -UmlU a (R) . (3) 

This equation is a consequence of the (exact) Heisenberg equations of motion. Now we 
specialize to nonrelativistic nuclear physics by limiting ourselves to that subspace of the 
complete Hilbert space that contains at most two (real) nucleons interacting by exchanging 
one (virtual) meson at a time. There the total Hamiltonian of the nucleus H is the sum of 
the kinetic and potential energies H = T + V of the nucleons, and the total axial current 
Jg(R) consists of one- and two-nucleon parts. We assume that the axial charge density p% [] 
is well approximated by its one-nucleon part | pi x _ h . This assumption agrees with - indeed it 
follows from - our fundamental assumption of nonrelativistic nuclear dynamics. This means 
that all (relativistic) operators are expanded in powers of 1/M; this provides a convenient 
book-keeping device in what follows. As shall be shown below, the (partial) continuity 
equation (|3|) strongly constrains (the longitudinal part of) the nuclear axial current, and in 
particular its two-nucleon, or meson-exchange part. 

We may break up the axial current conservation equation into one- and two-body parts 
without loss of generality. The divergence of the complete one-body current equals — % times 
the commutator of the kinetic energy T and the one-body axial charge density 

V • Jg(l - body) — —i[T, p 5 (l - body)] - / ff m=n a (l - body) , (4) 



2 In the following we shall drop the "index variable" R in the current and charge operators, except 
when necessary to avoid confusion. 

3 The subleading (relativistic) correction to the axial charge density operator contains two-body 
terms that modify the following analysis somewhat, but cannot change its main conclusion, due to 
an intrinsically different tensor structure of the one- and two-body operators. 
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is of 0(M~ 2 ), i.e., zero to leading order in 1/M, due to similar momentum dependencies 
of the kinetic energy T and the axial charge density — body) operators, as well as to 
the absence of non-diagonal isospin operators from T. Therefore, the test of conservation 
of the complete nuclear axial current is whether or not the potential V commutes with the 
one-body axial charge density. It turns out that, due to the momentum operator inside 
of pg(l — body), only a completely trivial, viz. a spatially everywhere constant potential 
commutes with the axial charge. In nuclear physics, therefore, one always needs a two-body 
axial current axial current Jg(2 — body) = J2f<k J5,(jfc)(2 — body) to compensate for the 
temporal change of the axial charge density. 

To show this formally we note that in general there are two possible sources of the non- 
commutativity of meson exchange potential V and the one-body axial charge density p%: (i) 
non-commuting isospin factors; and (ii) non-commuting spin-spatial factors, as can be seen 
from the identity 



V ■ J£(2 - body) = -i [V, p%] - hmiU a {2 - body) 



A A 



T, 



(0 



2M 



T 



2EE 

i=l j<k . 

• A A 

-^EE 

i=i j<k 

^(2 -body) 



(0 



{V(jk), (i) ■ V(i), 8(R - r(j) 
V(jk), {<T(i) • V(i),5(R- r W) 



}} 



(5) 



where the two-body potential V(jk) — I(jk)V{jk) is a product of its isospin and spin-spatial 
VufA parts. Since 



Vy fc ), {V (i) , 5(R - r w )}] oc 6(R - r w )^ [v (i) , V ijk) ] + 



(6) 



and the isospin anticommutator , I^k) } 7^ does not vanish for at least one value of 
a = 1, 2, 3, we see that the next-to-the-last line in Eq. (||) does not vanish, and thence the 
commutator of the potential and the axial charge density does not vanish for any potential Vjk, 
except the trivial one, i.e., a constant: Vjk = const. 0. In other words, axial meson exchange 
currents (MEC) J^(2 — body) are always necessary, so long as two nucleons interact and 
their dynamics can be described by Quantum Mechanics []. The same conclusions, of course, 
hold in nonrelativistic (NR) quark models, i.e., axial two-quark currents are always necessary 



4 This argument does not furnish a formal proof so long as all possible spin and isospin operators 
have not been examined for accidental vanishing of their anticommutators with cr^ and r^, 
respectively. Our proposition has been confirmed in all cases that we have studied so far, which 
cases constitute some of the most important parts of the nuclear NN potential. 

5 This does not mean, however, that we know how to construct axial MECs that obey PCAC 
for arbitrary NN potentials. Presently we know how to calculate only those axial MECs that are 
related to potentials that are based on chiral Lagrangian models. Methods ordinarily used for 
electromagnetic (EM) MECs do not seem to work here, cf. Ref. [Q]. 
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in NR quark models, as well. It is the compelling nature of this argument that makes it 
distinct from previous arguments along similar lines . 

This is perhaps a somewhat surprising result in view of the fact that the EM current 
conservation does not require MECs for many parts of the nuclear potential, e.g. for the 
contributions from the exchange of neutral mesons. Hence it will be our task to construct 
axial MECs associated with the exchange of the most important, i.e., the lightest, mesons. 
The 7r-exchange axial current has been known for some time , so we shall not repeat its 
derivation here. In this paper we shall concentrate on the isoscalar scalar (a) and vector (u) 
mesons. The isovector vector (p) meson cannot be introduced in a model-independent way, 
so we leave it for another occasion. 

To reveal the necessity of consistency between the nuclear wave functions and the nuclear 
pion production/absorption amplitude we must remember that the operator equation ([3]) 
describing PCAC is just a shorthand for the same statement about all axial current nuclear 
matrix elements {J^ 5 )fi- In momentum space the Heisenberg equations of motion lead to 

g*W/i = q-<J»>/i 

= q-(J!>/i-(P, pt))ti 

only if the initial and final states 1^/,$) are solutions to the nuclear Schrodinger equation 
H\^ n ) = E n \^ n ) (here q = E f - We used the abbreviation (A) fi = f \A\%). 

Thus, the nuclear pion absorption/emission amplitude (T*)fi must be evaluated using a pion 
absorption/emission operator T° that matches the nuclear dynamics leading to the nuclear 
wave functions \^f n ), if the result is to agree with PCAC. In the last line of Eq. (|7|) we used 
q 2 = qQ—c[ 2 , which is true for elastic scattering in the Breit frame, i.e., when q = Ef—Ei = 0. 

The idea of consistency of nuclear wave functions and MECs is not a new one: it has 
been a part of the electronuclear physics "folklore" for some time. The extension of this 
idea to the axial currents and pion production amplitudes seems to be less well known. In 
the following we shall apply this idea to nuclear axial currents in several variations of the 
sigma model with uj mesons. 

III. AXIAL CURRENTS IN SIGMA MODELS 

We shall use two of the sigma models already developed in Ref. [ffl]. As is well known, 
the pseudoscalar and the pseudovector, or gradient, ttNN couplings reduce to the same 
one-pion-exchange potential (OPEP) to leading order in nonrelativistic (NR) expansion. At 
first sight one might think that this implies identical axial MECs in the linear and nonlinear 
sigma models. This is not so because the linear sigma model OBEP includes a sigma- 
exchange potential as well, whereas the nonlinear one does not. Moreover, the difference 
between the "heterotic sigma model" and the other two sigma models persists, although in 
an unusual way: its (spatial) axial current is renormalized by a factor g^, though the axial 
charge is not. In the following we look at each model separately. Subsequently we add the 
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isoscalar-vector meson u, which plays an important role by providing short-range repulsion 
and is chirally invariant by itself. That brings about new terms into the axial current and 
the pion production operator in a way that is consistent with PC AC. 



A. The one-body current 

The one-body part of the nuclear axial current is the sum over all nucleons of the direct 
axial current and the the pion pole term, see Fig. 3(a) in Ref. Ufl. In configuration space 
this takes its usual non-relativistic form 



Jf (1 - body) =9Aj2f [* {{) - V R [ j^Z^ JJ<H R - r ») • W 
The axial charge density 

pg(l - body) = -ig A g ^ {cr {i] ■ V W , 5(R - r (0 )} (9) 

on the other hand is taken without the pion-pole contribution. This simplifying assumption 
is justified ex post facto by the fact that it does not prevent a successful construction of a 
consistent approximation. This does not mean that one cannot find, or should not search 
for approximations that avoid this assumption. 

The axial coupling constant g A is either unity, as in the linear sigma model, or its 
measured value 1.26 in the spatial part of the axial current Eq. @, and unity in the axial 
charge Eq. (^j of the heterotic sigma model of Ref. []]]. [For more on this, see Appendix 
[A].] We also neglect all nucleon electroweak form factors. This should be adequate for the 
purpose of describing the static properties of the nucleus. The nuclear matrix element of 
the direct one-body axial current is depicted in Fig. |l|(a), whereas the pion-pole part can 
be constructed by attaching the axial current "wavy line" to the external pion in Fig. |](a). 

a. Linear sigma model The i t/l -nucleon current in the linear sigma model, i.e., with a 
unit ("normalized") nucleon axial coupling, in momentum space reads 



T (i) 



J6, w (P<>Pi) = -ir 



(*) _ q ~~2^ 2 



(10) 



where q = p.- — Pi, and satisfies the nonrelativistic (static) version of the single-fermion axial 
Ward-Takahashi identity (WT id.) 



V-( P ;, Pj ;l-body) 



^ q 2 — m 2 



"- z (q^rlf) r - (p: ' Pj;1 - body) ' 

which follows from the Goldberger-Treiman (GT) relation M = gofw We use the same 
symbol for operators in configuration and momentum space. The second line on the right 
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hand-side (r.h.s.) of Eq. (|TT|) is the single-nucleon pion absorption operator multiplied by 
the divergence of the axial current factor /^m 2 and the static pion propagator (q 2 + rn^)~ l . 
We see that the pion absorption operator arises naturally from the divergence of the axial 
current. Both in this and in the heterotic sigma model the one-body axial charge operator 
Eq. (P) does not contain the factor g^, i.e., in momentum space it reads 



P5,(i)(Pi>Pi) = 



2M 



(12) 



for proof see Appendix [A]. 

b. Heterotic sigma model In the heterotic sigma model the GT relation is modified to 
QaM = g^NNln, where gA = 1-26 and the one-body current becomes 



(0 



h,(i){PvPi) =9a^ 



(0 ~~ q ^2~i 2 

\ q^ + m~ y 



which satisfies the single-nucleon axial Ward-Takahashi identity (WT id.) 



(13) 



q ■ J 5,(i) (Pi> Pi) - A I 2 + ' r m 2 j dvNNTf-) ( 1 



2M 



r»(p;, Pj ;l-body) 



,q 2 + <, 



r»(p;, Pi ;l-body) 



(14) 



since [ T, pg(l — body)] = + (9(M~ 2 ). The spatial parts of one-nucleon axial currents in 
the two models are almost identical in the NR limit, the only difference being the overall 
factor gA- As stated above, the axial charges in the two models are identical, (see Appendix 

0)- 



B. Two-nucleon axial current 



Having shown that in order to have a partially conserved axial current in a non-relativistic 
nuclear model, one must have two-nucleon axial currents, we turn towards constructing such 
MECs. The two-body currents Js(2 — b) appropriate to the meson exchange two-nucleon 
potentials will be constructed and we shall show that they lead to PCAC. We construct these 
"meson exchange currents" by non-relativistic reduction of covariant Feynman amplitudes 
in specific chiral models. 

Due to the presence of gradient operators it will be to our advantage to work in the 
momentum space. Definition of the Fourier transform of two-body currents into momentum 
space can be found in Ref . , among others. The current conservation relation in momentum 
space can now be written as 

q ■ Jl(2 - body) = [V(2 - body), - body)] - i [ K(2 - body) . (15) 
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We repeat that this equation follows from only two assumptions: (i) PCAC, and (ii) quantum 
mechanics. As stated above, the only "sure-fire" way of constructing an axial MEC that 
satisfies PCAC that we know of is to start from a relativistic chiral Lagrangian model. We 
shall use the two variations of the linear sigma model already utilized in Ref. [|IJ and the 
simplest uNN interaction Lagrangian that preserves chiral symmetry. 



1. Linear sigma model 



To construct the partially conserved nonrelativistic axial two-nucleon current in this 
model we start from the corresponding covariant amplitude. A nonrelativistic expansion in 
powers of 1/M leads to three terms of 0(1/ M): (i) one due to the meson-in-flight diagrams, 
Fig. (0) plus their pion-pole counterparts, Fig. (|J), (the covariant current that defines this 
amplitude can be found in Eqs. (33), (34) of Ref. Jl|): 



A?(k!,k 2 ,q) 



9l a 

lT(- 



2M (1) 



k x - k 2 + q 



f-KQo 



X 



<T(i) ■ ki 



(ki 



v q 2 + m 2 J 
+ (1 <-> 2) 



(16) 



m 2 ) (k 2 + ml) 

J - ml. Three-momentum conservation reads 



where k; = p; - p., i = 1, 2; and Ug a ^ = m^ 
\ti + k 2 + q = 0. This MEC alone is not sufficient for PCAC, as can be seen from the 
corresponding divergence, which reads 



q • A?(k l9 k 2 , q) = jLfa [kg - k 2 + (ml - mj) - Um\ 



9o 



v q 2 + ml 



<T ( i) • ki 



(k 2 + ml) (kf + ml ; 

.a a (1) ' k l 



9l 



2M 



+ 



+ 



+ 



A. 

2M 

2M 
2M 



(1) k 2 + m, 

a ■ k l 

kf + ml 



+ (1 <-> 2) 
- (1 «-> 2 



+ (1 - 2) 



+ O(Uml) 



'«kl 



m 2 



2) 



where the one-meson-exchange NN potential in the linear sigma model reads 

9l 

p 2 + ml 



(17) 



S 



+ T (1) 




(18) 



Note that the cr-exchange potential is of 0(1/M° = 1), whereas the 7r-exchange potential is 
of 0(1 /M 2 ). Consequently, the commutator of the potential and the one-body axial charge 
also fall into two distinct orders in 1/M. This fact allows an apparently clean and simple 
separation of the cr-exchange current effects from the 7r-exchange ones. In this paper we are 
primarily interested in the cr-exchange currents, so we leave the 7r-induced ones aside, as 
they have been extensively studied in the literature ||||. 
The commutator on the right-hand side of Eq. (|1| 



is 



2-b, 



[V a , pl\ + 0(1/M 3 



2M 



T, 



'(1) 



•k. 



(1) k| + ml 



+ (1 



+ 0(1/M 3 



(19) 



We see that the cr-exchange leads to terms of 0(1/M), whereas the 7r-exchange leads to terms 
of 0(1/M 3 ) in the commutator. Comparing Eq. flT9"| ) with the divergence of the axial two- 
body current Eq. ([[7]) we see that we are rather far from having the sigma-exchange potential 
commutator on the right-hand side. There are four terms left over: two are proportional 
to the 7r propagator thus indicating perhaps a relationship to the 7r-exchange current, but 
also being of 0(1/ M), i.e., two orders in 1/M lower than the lowest expected 7r-exchange 
current contribution! The other two are proportional to the cr propagator. In other words, 
our expectations do not square well with these initial results. So, the question is: whence do 
these "extra" terms come from and is there something that might compensate for them? The 
answer is that among the apparently "higher-order" terms in the axial 7r-exchange currents 
there are two such ones. 

(i) One is the 7r-exchange axial current, Fig. f|(b), 



An^k^q) 



9o 



/ttQ 



2M 2 \q 2 + ml 



(1) k[T 



mi 



(20) 



that is actually one order in 1/M lower than naively expected, due to the validity of the 
Goldberger-Treiman relation f n g = M in the linear sigma model. Thus 



An^k^q) 



9l 



q 



2M V q 2 



mi 



r a 0-(l)-kl 

f(1) kTT 



mi 



2) 



(21) 



This current is due to an effective two-pion-nucleon vertex, which is really a time-ordered 
Z-graph arising in the NR reduction of the nucleon Feynman propagator (the so-called pair 
current"), and is not an "elementary" interaction term in the Lagrangian. The divergence 
of this axial MEC is 



q- A&(ki,k 2 ,q) 



9o 



2M 



a q-(i) • k i 

{1) kf + m 2 



+ (1 



+ O(Uml) . 



(22) 



^ii) Similarly, there is a cr-exchange "Z-graph", Fig. ^|, as well 
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q 



2M I q 2 + m 2 



a o-(i) ■ q 

(1) k 2 + m 2 



+ (1 <->2) 



(23) 



Its divergence is 



q- A? n (ki,k 2 ,q) 



2M 



(1) (k|.+ 



^ + (1-2) 



+ O(Uml) 



(24) 



which is exactly what we need to bring the continuity equation (|i~7|) in the linear sigma 
model into the canonical form flTo]). 

Thus, the complete nonrelativistic axial MEC to 0(1 /M) in the linear sigma model is 

J5_ b (ki, k 2 , q) = Af (k 1? k 2 , q) + Af^kj, k 2 , q) + A^kx, k 2 , q) 

CTTVTT 




2M {1) 



ki 



q 



q 2 + m 2 



x 



(i) 



ki 



(k 2 + m 2 ) (k 2 + m 2 

( q 



q 2 + m 2 



a (l) • ki _^ • q 

k? + m 2 k| + m 2 



+ (1 <-> 2) . 



(25) 



Thus, Eq. ( |25| ) uniquely fixes the pion production/absorption operator in the one- 
boson-exchange potential approximation to the linear sigma model. This operator is to 
be used together with the 0(1/M 3 ) 7r-exchange operators and nuclear wave functions that 
are solutions to the Schrodinger equation with the above one pion + sigma exchange NN 
potential. 

The resulting MEC operator is perhaps somewhat unexpected: certainly the first 
( "meson-in-flight" ) term is not a surprise, but the presence of the second ( "pion seagull" ) 
term might seem a little odd at first: One would have been hard pressed to correctly guess 
the second term in the MEC without the benefit of guidance by the linear sigma model. 
Thus we have expanded the nuclear interaction to include one-sigma-exchange NN potential 
(beside the OPEP) in a manner that is consistent with PCAC, but still with qa = 1. Next, 
we shall relax that assumption. Once again, we resort to a chiral Lagrangian model for 
guidance. 



2. Heterotic sigma model 

The heterotic sigma model differs from the linear one by the presence of gA in the GT 
relation, g n NNfn — 9aM, and in the spatial part of the axial current (||D, but not in the axial 
charge (0). [For a proof of this statement, see Appendix [A].] This variation induces some 
curious changes in the axial two-nucleon currents. To begin with, there are two new types 
of relativistic Born approximation Feynman diagrams, depicted in Figs. [TJ(b), f|(b), and Fig. 
|6]. [The complete relativistic current consists of the sum of Eqs. (41), (47) and (51), and 
its divergence in the sum of Eqs. (42), (48) and (52) in Ref. []I]].] By evaluating its matrix 
element between on-shell single-nucleon states and making the nonrelativistic expansion in 
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powers of 1/M, one finds a host of 7r-exchange currents of 0(1/M 3 ), and two different o- 
exchange currents of 0(1/M). One of them is the familiar meson-in-flight diagram (+ its 
exchange), Fig. |2], but rescaled by a factor of g^' 



Af(k 1 ,k a ,q;/ l )=p x ^7^ ) 



ki - k 2 + q 



f*9o 



<T(i) ■ ki 



(k 2 + m 2 )(k 2 + m 2 
for which the current divergence in momentum space reads 



q 2 + m 2 
+ (1 <-> 2) , 



q-Af(k 1| k 2l q;/i)=^^7f 1) 



Ug 



X 



(Tii) ■ ki 



(k| + m 2 ) (k 2 + m 2 ) 



q 2 + m 2 / 
+ (1 <-> 2) 



<? 2 


T (l) 


(T(i) • k x 


+ (!<- 


-2) 


2M 


kf + m 2 


So 2 
2M 


T (l) 


o-(i) • q 

k| + m 2 


+ (1- 


-2) 



+ 9A- 

+ 9a: 
+ 0(f n ml) . 

The one-meson-exchange NN potential in the heterotic sigma model reads 

V^-b(p) =K(p) + K(p) 

9l 



p2 _|_ m 2 



_ _ /<T(l) • p\ /CT (2 ) • p\ 



91nn 
P 2 + 



9l 



(26) 



(27) 



(28) 



and c/nNN = <?vi<7o> due to the new GT relation. 

Similarly, the pion "Z-graph", Fig. |l|(b), contribution is also renormalized by factor g&: 



An(ki,k 2 ,q; h) = -g A 



9l 



qA 



-9a 



2M 2 \q 2 + ml t 

q \ 



9l 



2M \q 2 + ml / 
The divergence of this axial MEC is 



r «k 2 + m 2 +(1 " 2) 

^)S^ + (l-2) 
w kf + m 2 



(29) 



q- A5(ki,k 2 ,q; h) = -g A 



9t 
2M 



a <T (1) ■ ki 

«k 2 + m 2 +(1 ^ 2) 



+ 0C/X 



(30) 
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Finally, in the heterotic sigma model in addition to the cr-exchange "Z-graph", Fig. 
contribution 



A^ II (k 1 ,k 2 ,q;/ i ) = -||:r { ' 1 1) 



q 



\ °"(i)-q 

q 2 + m 2 / k 2 , + m 2 

there is also an elementary A^aNN axial current vertex, Fig. | 
MEC, Fig. |, 



+ (1 



(31) 



and a naNN vertex-induced 



ArV(ki,p 2 ,q; h) = -/ w 



'9a -I 



9l 



'(i) 



U /2M(k 2 + m 2 ) 



2<7(i) 



q 



cr (1) • (k 2 - q) 



q 2 + m 2 j 

The divergence of the sum of these last two axial MECs is 



+ (1 <-> 2) . 



(32) 



q- A in+iv( k i. k 2>q; h) 



9l 



T, 



(1) 



(g A - l)o- ( i) • (k 2 + q) + o-(i) • q 



2M kg + m 2 

+ (1 <-> 2) + 0(Um 2 n ) , (33) 

which is exactly what we need to complete the continuity equation fllBl). 

The sum of these three axial MECs together with the one-body (impulse approximation) 
terms plus the 0(1/M 2 ) pion-exchange currents which we consistently suppressed in this 
paper, constitute the complete, PCAC-obeying axial current: 

a a _,_ a a _| a a _i_A a 
^-I ' A ll i ^-III i A IV 

ki - k 2 + q 



J 5,2-b( k i ; k 2 ; q; h ) 



9 T a /„ 



q^ + 



x 



r (i) 



ki 



(k 2 + ml) (k 2 + ml) 



q 



q 2 + m 2 t 
' 9a -I ^ 



9a 



& (l) • ki | <T(i) • q 
kf + m 2 



+ 



k 2 + m 2 



2<r (1) + 



q 



q 2 + m 2 7 



(i) 



■ ( k 2 - q) 



+ (1 



(34) 



Most of this current is new: specifically, terms (II, III, IV) have not been considered before, 
only the first term (I) having been derived in Refs. [|K],|r]]]. This current satisfies PCAC 

q ■ J^ b (k 1; k 2 , q; h) = [V a , p a 5 ] + 0(f n ml) . (35) 

This equation uniquely fixes the pion production/ absorption operator in the OBEP approx- 
imation to the nuclear dynamics in the heterotic sigma model as the u O(f n mD" term, to 
be evaluated in the next section. We see that, by a curious turn of events, the divergence 
of the complete axial current (|35|) equals the divergence found in the linear sigma model 
(p7|), modulo different 0(f n ml) terms, of course. In other words, we have exactly the same 
commutator of the one-sigma-exchange potential (OSEP) and the axial charge in these two 
sigma models, despite manifest differences between their axial currents. This fact is a con- 
sequence of identical axial charges and OSEP in the two models, the former being a subtle 
effect explained in [A]. Finally we turn to the u exchange. 
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3. Omega meson exchange 



One may include the u;-meson-exchange potential 



Vu, 



u: 



p 2 + m 2 



(36) 



into the nuclear Hamiltonian H at no peril to chiral symmetry because it is an isoscalar vector 
field which is chirally invariant by itself in both the linear and the nonlinear realization. As 
is well known, the main benefit of including this term into the nucleon-nucleon potential is 
that it provides short range repulsion that is otherwise absent from the sigma models. The 
nuclear potential with ir, a, u exchange has sufficient attraction to bind the deuteron and 
enough repulsion to keep it weakly bound. The associated axial MEC is, Fig. [7|, 



J L ( k i, k 2,q) = -f*g* 



nn: 



91 



q 



X 



2M 2 \q 2 + ml 
C(i) • k 2 



(1) kg+m£ 



+ 1-2 



(37) 



This completes the construction of the PCAC-constrained axial current in models with 
OBEP based on the n, cr, and uo mesons. The p meson was deliberately omitted since its 
contribution is (highly) model dependent. 



IV. PION PRODUCTION OPERATORS 

The one-body pion production operator is well known, as well as the two-body one asso- 
ciated with pion exchange. We shall therefore concentrate on the MECs that are associated 
with other meson (a and u) exchanges, as specified by the PCAC constraint Eq. (|f5|). 



A. Linear sigma model 

It follows from Eq. flT5| ) and the linear sigma model axial current Eq. 

r« 2 _ b (k x , k 2 , q) = I* (k l5 k 2 , q) + r« n (k 1; k 2 , q) + r« m (k 1; k 2 , q) 



that 



■ So 3 



2M 2 W 

(T7T7T 



Q"(2) • k 2 <T( 2 ) ■ q 
k 2 + ml kf + ml 



9l 



2M' (2) (k? + ml ) (k| + ml 



(2) 



+ (1 - 2) 



(38) 



as the complete linear sigma model pion production operator to this order in 1/M. This 
result corresponds to Figs. f|(b), || |^. Although it has long been known that the corre- 
sponding covariant amplitude is chirally symmetric, we are not aware of anyone having used 
the MEC ( |3*8"D in nuclear physics. 
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B. Heterotic sigma model 



As stated before, the main difference between the linear and the heterotic sigma models is 
the axial coupling constant g^, which is induced by the new derivative-coupled interactions. 
The latter, in turn, renormalize one old graph, Fig. [|, and creates two new elementary 
diagrams: Fig. |6], and Fig. ||(b) [J Thus one finds 



2 _ h (k 1; k 2 , q; h) = T a nl (k 1} k 2 , q; h) + T a nU (k 1; k 2 , q; h) 
+ I'",,, lk,.k.,.( ' ' 1 



,ri k 2 , q; ri) + 1 v n iki, K2, q; n 

: m (k^k^q; h)+T a nW (k^k^q; 

• ffo « [ o-(2) • k 2 (7(2) • q 
" 2M 2 (2) [ 5A k| + m 2 ^ kf + m 2 , 

+(^-i) ^ 2 ^; q) 

kf + m% 

-jn n , 9 T a 
igcnrirgA >> 



h) 



9o a ° (2) ■ k 2 

2M (2) (kf + ml ) (kl + 



+ (1-2 



(39) 



This is the complete two-body pion-production operator in the heterotic sigma model, 
though with all 0(1/ M 3 ) pion-exchange two-body operators excluded. This operator is 
only to be used with nuclear wave functions that are solutions to the Schrodinger equation 
with the aforementioned one pion + sigma exchange NN potential, Eq. 



C. Omega meson exchange 

The omega-exchange 7r-production operator, Fig. 0, to this order in 1/M is 

KM , q) ^ ls f^^^l + {1 ^ 2) ). (40) 

Note the utter functional dissimilarity between this and the two a-MECs (^8j), and (j39l), 
despite the fact that both are due to terms in the NN potential that are indistinguishable to 
this order in nonrelativistic expansion, viz. the isoscalar one-vector- and the scalar-meson 
exchange potentials in Eqs. (0), and fl28| ) respectively. Specifying the Lorentz structure 
of the NN potential distinguishes the vector- from the scalar-exchange induced terms, but 
even then the ambiguity between the two kinds of scalar MECs is survives. This fact brings 
into focus the intrinsic and apparently intransigent ambiguities associated with attempts to 
construct consistent axial MECs starting from the NN potential [^] . At this stage the only 



Scalar meson exchange currents have been introduced into the analysis of pp — > irpp on an ad 
hoc basis [QJ5). Those papers are different from ours in that they do not include the pion-sigma- 
exchange graph, nor was there any concern shown for the consistency between the nuclear wave 
functions and the pion production operators. 
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formalism that allows systematic construction of consistent is based on relativistic chiral 
Lagrangians. 

We believe the preceding results to be important for two reasons: (i) They provide several 
explicit examples of how PCAC taken as an underlying principle leads to classification 
and construction of admissible (PCAC-consistent) approximations to nuclear axial current 
matrix elements in chiral models, (ii) They specify PCAC-consistent, apparently novel pion 
production operators in several OBEP nuclear models. Besides the pp — > n pp reaction, 
|13fl , these results ought to also be directly applicable to studies of the pp — > 7r + d reaction 

V. SUMMARY AND CONCLUSIONS 

In summary, in this paper we have used PCAC to relate the nuclear axial current and 
pion production operators to each other, and to the two-nucleon potential, in nuclear theories 
based on the nonrelativistic Schrodinger equation. We focused in particular on the axial- 
and pion-production MECs related to the exchange of the lightest isoscalar mesons with 
J p = + , 1~, i.e., to a, uj mesons. [Pions have been treated elsewhere |8|||, and the p 
meson contributions are highly model-dependent.] We constructed axial currents and pion 
production operators in two variations of the linear sigma model, with = 1, or g& = 1.26, 
with or without cj-exchange and showed explicitly that they satisfy the PCAC constraints. 
All of these results are new, to our knowledge. 

In the process we also made several assumptions and approximations: 

(i) we neglected the pion-pole term in the axial charge. This is justifiable since no need 
arose for it in our analysis. 

(ii) we neglected all retardation effects, as well as the recoil MECs. 

(iii) we neglected all isovector- vector and/or axial vector meson exchanges. 

(iv) we did not include nucleon or meson form factors, either electro-weak or strong. 
Another objection, perhaps of a more theoretical nature, can be raised against the present 

calculation: to talk about meson production in a NR potential theory (where such meson 
d.o.f. have been "integrated out") seems self-contradictory. We have side-stepped this 
problem by defining the pion production operator as the O^f^m^) term in the divergence 
of the axial current Eq. ([!]), in analogy with the relativistic PCAC result. This issue can 
be addressed more deeply within the Fukuda-Sawada-Taketani-Okubo-Nishijima (FSTON) 
method for constructing effective nuclear theories [T5|], and we are currently working on it. 



Partial conservation of the axial current ought to be an important criterion in the con- 
struction of both nuclear two-body, or meson-exchange axial current operators, and of nu- 
clear two-body pion-production operators. The latter arrise from two sources: (i) "elemen- 
tary" meson-nucleon vertices present in the chiral Lagrangian; and (ii) "effective" meson- 
nucleon vertices due to the time-ordered Z-graphs, also known as pair currents. They have 
not, to our knowledge, been examined from the present viewpoint, with the exception of 



7 This includes the nuclear chiral perturbation theory {N\PT), which is in the present OBE 
approximation just the one-pion-exchange term that was previously treated in Refs. |||f| 
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Refs. HH on the one-pion-exchange axial currents and the early work by Blin-Stoyle and 
Tint 0] on nuclear pion production. Our study is only the first step in joining these two 
ideas and extending them to include light isoscalar mesons (<j,uj). 

This work is based on two fundamental assumptions: (i) quantum mechanics; and (i) 
PCAC. Although our examples were drawn from two specific sigma models, they exemplify a 
far more general relation between the nuclear Hamiltonian and the nuclear pion production 
operator. Indeed such a relation must hold in any calculation based on the two aforemen- 
tioned principles, and in nuclear chiral perturbation theory in particular. This relationship 
is almost completely unexplored in the last mentioned setting, a situation that must be 
remedied. 

Note added in proof. We have learnt of the related work by S. M. Ananyan, Phys. Rev. 
C 57, 2669 (1998) only after the present paper had been submitted. 
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APPENDIX A: AXIAL CHARGE DENSITY IN THE HETEROTIC SIGMA 

MODEL 



The heterotic sigma model is a chirally symmetric field theoretic model that leads to 
an axial current with arbitrary Qa(^ 1) 5 and a mixture of pseudoscalar and pseudovector 
pion-nucleon couplings ffl. The Lagrangian density of this model is given by 



C = ijipifj - gotfj [a + i l5 7z ■ r] i, + l - (d,4>f - V{<t> 2 ) 



+ 



' 9a ~ 1 

v fir 



(Al) 



where </> = (a, 7r) is a column vector and V is the same potential as in the linear sigma 
model. The (partially conserved) axial-vector Noether current in this model reads 



J IT 



+ 



+ o- [ V>7m75^ ) +a y^t^^ x 7r 



(A2) 



After shifting the sigma field this can be written as 
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TO 



+ 



9a (^1^15 ^ ) + U^ir + (sd^n - 7T<9 M s)° 
' 9a - 1^ 



J TV 



+ s {2f n + s) ^7m75tt^ +(fn + s)[ ipi^-i) x 



7T 



(A3) 



Thus we see that the nucleoli axial current has acquired a new coupling constant: 1, 
which was the purpose of this model. 

The (new) derivative coupling terms in the Lagrangian ( |A1| ) modify the canonical mo- 
menta as follows 



n CT = & - j 

IE = 7T a + 



, ft J 

( 9A-1 

v /. 2 , 



/t T ' 77 / 



, + T x 77 , \ ,+ ' , 



(A4) 
(A5) 



We see that the axial charge retains its linear sigma model form when written out in terms 
of canonical fields and their associated momenta: 



pl = JS 5 = ^75 - (^ a n CT - an*) . 



(A6) 



Hence we see that the axial charge carried by the nucleon is unchanged as compared with 
the one in the linear sigma model, i.e., we have 9a — 1 here, which was to be proven. 



17 



REFERENCES 



[1] V. Dmitrasinovic, Phys. Rev. C 54, 3247 (1996). 

[2] R.J. Blin-Stoyle and M. Tint, Phys. Rev. 160, 803 (1967). 

[3] D. Koltun and A. Reitan, Phys. Rev. 141, 1413 (1966). 

[4] C.J. Horowitz, Phys. Rev. C 48, 2920 (1993); C.J. Horowitz, H. O. Meyer and D. K. 
Griegel, Phys. Rev. C 49, 1337 (1994); see also J.A. Niskanen, Phys. Rev. C 53, 526 
(1996). 

[5] T.-S. H Lee and D. O. Riska, Phys. Rev. Lett. 70, 2237 (1993). 

[6] see e.g. R.J. Blin-Stoyle, p. 3 - 66 in vol. I, and J.-I. Fujita, and M. Ichimura, p. 625 - 

650, in vol. II, eds. M. Rho and D. Wilkinson, Mesons in Nuclei, 
[7] Y. Futami, J.-I. Fujita, and H. Ohtsubo, Suppl. Prog. Theor. Phys. 60, 83 (1976). 
[8] K. Ohta, Nucl. Phys. A 314, 333 (1979). 

[9] J. Adam jr., Ch. Hajduk, H. Henning, P. U. Sauer and E. Truhlik, Nucl. Phys. A 531, 
623 (1991); and I. S. Towner, Nucl. Phys. A 542, 631 (1992). 
[10] K. Tsushima and D.O. Riska, Phys. Lett. B 333, 17 (1994). 
[11] S. Schramm and C.J. Horowitz, Phys. Lett. B 335, 279 (1994). 
[12] P. G. Blunden and D.O. Riska, Nucl. Phys. A 536, 697 (1992). 
[13] T. Sato, T.-S. H. Lee, F. Myhrer, and K. Kubodera, Phys. Rev. C 56, 1246 (1997). 
[14] L. Canton, A. Davini, and P.J. Dortmans, preprint nucl/th-9801028. 
[15] M. Kobayashi, T. Sato, and H. Ohstubo, Prog. Theor. Phys. 98, 927 (1997). 



18 



FIGURES 



FIG. 1. Effective nonrelativistic Feynman diagrams contributing to the one- (a) and the 
two-body axial current nuclear matrix element (b). Each graph consists of a "direct" term and a 
"pion pole" term. We display only the direct terms here. In the nonlinear sigma model there are 
only graphs 1(a) and 1(b). The dashed line denotes a pion, the solid one a nucleon, the wavy line 
is the external axial current (source). 

FIG. 2. The "meson- in-flight" graph contributing to the axial current in the linear and the 
heterotic sigma models. The zig-zag line denotes a sigma meson. One must keep all four graphs 
in Figs. l(a),(b) and 2, 3 in both the linear and the heterotic sigma model. 

FIG. 3. Sigma-meson-exchange current contributing to the nuclear pion production matrix 
element. It consists of a "pair", or Z-graph time-ordered contribution in the linear sigma model to 
which the elementary iraNN vertex is added in the heterotic model. 

FIG. 4. Effective nonrelativistic Feynman diagrams contributing to the one- (a) and the 
two-body pion production nuclear matrix element (b). In the nonlinear sigma model there are only 
graphs 1(a) and 1(b). The pion-exchange current (b) consists of a "pair", or Z-graph time-ordered 
contribution in the linear sigma model to which the elementary ttttNN vertex is added in the 
heterotic model. The dashed line denotes a pion, the solid one a nucleon. 

FIG. 5. The "meson-in-flight" graph contributing to the linear and the heterotic sigma models. 
The zig-zag line denotes a sigma meson. One must keep all four graphs in Figs. l(a),(b) and 2, 3 
in both the linear and the heterotic sigma model. 

FIG. 6. Sigma-meson-exchange current contributing to the nuclear pion production matrix 
element. It consists of a "pair", or Z-graph time-ordered contribution in the linear sigma model to 
which the elementary iraNN vertex is added in the heterotic model. 

FIG. 7. The w-exchange current contributing to the nuclear pion production matrix element. 
The curly line denotes an omega meson. This graph can be added to all sigma models without 
disturbing their chiral symmetry. 
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